Hidden magnetic order in CuNCN 



(N 

o 

(N 



(N 



C/3 



Alexander A. Tsirlin/^LJ Alexander Maisuradze,^ Jorg Siehelschmidt,^ Walter Sehnelle/ 
Peter Hohn.^ Ronald Zinke,'^ Johannes Richter,'^ and Helge Rosner^'IU 

^ Max Planck Institute for Chemical Physics of Solids, Nothnitzer Str. 40, 01187 Dresden, Germany 

^Laboratory for Muon Spin Spectroscopy, Paul Scherrer Institute, 5232 Villigen PSI, Switzerland 

^Institute for Theoretical Physics, University of Magdeburg, P.O. Box 4120, 39016 Magdeburg, Germany 

We report a comprehensive experimental and theoretical study of the quasi-one-dimensional quan- 
tum magnet CuNCN. Based on magnetization measurements above room temperature as well as 
muon spin rotation and electron spin resonance measurements, we unequivocally establish the local- 
ized Cu+'^-based magnetism and the magnetic transition around 70 K, both controversially discussed 
in the previous literature. Thermodynamic data conform to the uniform-spin-chain model with a 
nearest-neighbor intrachain coupling of about 2300 K, in remarkable agreement with the micro- 
scopic magnetic model based on density functional theory band-structure calculations. Using exact 
diagonalization and the coupled-cluster method, we derive a coUinear antiferromagnetic order with 
a strongly reduced ordered moment of about 0.4 ^ib, indicating strong quantum fluctuations inher- 
ent to this quasi-one-dimensional spin system. We re-analyze the available neutron-scattering data, 
and conclude that they are not sufficient to resolve or disprove the magnetic order in CuNCN. By 
contrast, spectroscopic techniques indeed show signatures of long-range magnetic order below 70 K, 
yet with a rather broad distribution of internal field probed by implanted muons. We contemplate 
the possible structural origin of this effect and emphasize peculiar features of the microstructure 
studied with synchrotron powder x-ray diffraction. 



PACS numbers: 75.50.Ec, 76.75.+i, 76.30.Fc, 75.10.Pq 
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I. INTRODUCTION 

While the majority of magnetic systems develops long- 
range order (LRO) of spins at low temperatures, the lack 
of LRO down to zero temperature is a less common sce- 
nario that implies strong quantum fluctuations and pe- 
culiar phenomena related to a specific type of spin lat- 
tice and/or exchange couplings. Low-dimensionality and 
frustration, two main prerequisites of quantum fluctua- 
tions, impede LRO in a quantum spin system by reduc- 
ing the nragnetic ordering temperature, the ordered mag- 
netic moment, and other features of the ordered statejii^ 
Although not necessarily preventing an ordered arrange- 
ment of spins at K, quantunr fluctuations render the 
LRO barely visible for experimental techniques that give 
response proportional to the size of magnetic moments 
(neutron scattering), transition entropy (specific heat), 
or other quantities related to the LRO state. For ex- 
ample, Sr2Cu03, a paradigmatic quasi-one-dimensional 
spin-i magnet, features the nearest- neighbor exchange 
coupling of J ~ 2600 K^^ yet undergoing LRO at as 
low as 5.4 K {Tn/J ~ 0.002) with the ordered mo- 
ment of 0.06 /is, compared to 1 /.Js in a classical spin- 
i system4. The LRO in Sr2Cu03 can only be observed 
with single-crystal neutron scattering^ or spectroscopic 
methods, such as muon spin rotation (/iSR)^!^ and nu- 
clear magnetic resonance (NMR))^ while powder neutron 
diffraction essentially fails to detect the LRO because of 
the vanishingly small ordered moment f2, 

The presence of a low-temperature LRO in some other 
low-dimensional spin systems remains a matter of con- 
tention. It is often declared that the lack of any visi- 
ble anomalies in thermodynamic properties^ as well as 



the lack of the magnetic neutron scattering in a powder 
experiment)^ are a good evidence for the absence of LRO. 
However, /tSR does find the magnetic transitions in many 
low-dimensional systems previously considered magneti- 
cally short-range-ordered in the experimentally accessible 
temperature range.— i^iii^ A reliable study of the LRO in 
a low-dimensional spin system requires a comprehensive 
experimental investigation combined with a sound micro- 
scopic analysis that provides details of the LRO, such as 
the magnetic structure and the anticipated ordered mo- 
ment. The microscopic input is essential for interpreting 
the experimental data, because the lack of the expected 
signal (e.g., in a powder neutron diffraction experiment) 
may either indicate the zero ordered moment (the ab- 
sence of LRO), or simply show that the LRO state is 
beyond the sensitivity threshold of the method. In the 
following, we apply a combined experimental and theo- 
retical approach to explore the ground state of CuNCN, 
a quasi-one-dimensional spin- 2 quantum magnet contro- 
versially discussed in the recent literature. 

CuNCN is a semiconducting compound containing 
spin-i Cu+^ cations and anionic carbodiimide [NCN]^~ 
groups. The layered crystal structure of CuNCN is 
formed by characteristic CUN4 plaquettes that share 
edges and form chains along the a direction. The nearly 
linear NCN groups connect the structural chains along c 
(Fig. [TJjii An early experimental study by Liu et alt^ 
reported the lack of LRO, as concluded from thermo- 
dynamic measurements and neutron powder diffraction. 
The disordered ground state was understood as a "highly 
correlated antiferromagnetic state at room temperature" 
caused by the frustrated triangular spin lattice formed 
by the couplings Ji and Ji,^ Later, the same group 
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FIG. 1. (Color online) Top panels: crystal structure of 
CuNCN, open circles denote Cu atoms centering CuN4 pla- 
quettes (shaded). Bottom panels: microscopic magnetic 
model proposed in Ref. [isl : open and filled circles at the sites 
of the spin lattice show up and down spins in the columnar 
AFM ground state (see Sec. lIVA|) . 



performed a polarized neutron scattering experiment^ 
that did not show any signatures of LRO, either. How- 
ever, this observation was now interpreted as a non- 
magnetic state of Cu"*"^. In a subsequent theoretical 
study, Tchougreeff and Dronskowskii^ recalled the sce- 
nario of the magnetic Cu"*"^ cations, and developed a 
mean-field theory for the resonating- valence-bond (RVB) 
ground state of an anisotropic J\ — Jf, triangular spin lat- 
tice (see the bottom left panel of Fig. [IJ. 

In a preceding study^i^ we proposed an independent 
and essentially different microscopic picture^ based 
on conventional density functional theory (DFT) band 
structure calculations that are well known as a highly ef- 
ficient tool for unraveling complex spin lattices in quan- 
tum magnets fi^ The leading antiferromagnetic (AFM) 
exchange J ~ 2500 K was found along the c direction 
via the NCN groups (Fig. [IJ. A weaker FM coupling 
J\ ~ —500 K connects the resulting spin chains along 
a, whereas a couple of marginal AFM couplings induce 
a weak frustration that should not preclude the LRO 
in CuNCN. In this paper, we provide experimental and 
further theoretical support for this microscopic scenario. 
Our data clearly indicate the presence of localized mag- 
netic moments on the Cu"*"^ atoms, and suggest the for- 
mation of static magnetic fields below 70 K. To under- 
stand the latter observation, we investigate the ground 
state of the proposed magnetic model, evaluate the Nccl 
temperature, and explain why the anticipated magnetic 
order has likely been overlooked in the previous stud- 
ies. Note that recently Zorko et alw^ presented another 
experimental work on CuNCN and proposed an "inho- 
mogeneous magnetic ground state" possibly related to 
the strong frustration on the triangular spin lattice. Al- 
though most of their experimental results match our find- 



ings, the interpretation is notably different. We discuss 
these differences in Sec.|V]of the present paper and in a 
separate Comment J^ 

The outline of our paper is as follows. We start with 
methodological aspects in Sec. [ill and proceed to the ex- 
perimental results in Sec. IIIII In Sec. IIVI we present a 
theoretical study of the microscopic magnetic model pro- 
posed in Ref. Il5|. and re-analyze the neutron-scattering 
data from Refs. [Ij and[l3l. We complete our work with 
a discussion and summary in Sec. |Vl 



II. METHODS 

Poly crystalline samples of CuNCN were prepared by 
a two-step procedure described in Ref. [III. We used 
CuCl2 • 2H2O and H2NCN as starting materials and 
Na2S03 as a reducing agent. H2NCN was stored in a 
fridge to avoid polymerization. CuNCN was obtained in 
the form of black powder. Repeated attempts of crys- 
tal growth were so far unsuccessful. Although labora- 
tory powder x-ray diffraction (XRD, see below) showed 
similar patterns for all samples, further studies (chemi- 
cal analysis, synchrotron XRD, magnetization measure- 
ments) revealed deviations in the chemical composition, 
microstructure, and magnetic behavior. 

Powder samples of CuNCN were characterized with 
laboratory XRD (Huber G670 Guinier camera, CuKd 
radiation, Ge monochromator, ImagePlate detector, 
26* = 3 — 100° angle range) and conventional analyti- 
cal techniques fi^ High-resolution XRD patterns for se- 
lected samples were measured at room temperature at 
the ID31 beamline of European Synchrotron Radiation 
Facility in Grenoble (wavelength A = 0.4 A, eight scin- 
tillation detectors preceded by Si (111) analyzer crystals, 
20 = 1 — 40° angle range). Samples were loaded into 
thin-walled borosilicate glass capillaries and spun during 
the data collection. 

Magnetic susceptibility was measured with a Quantum 
Design MPMS SQUID magnetometer in the temperature 
range 2 — 650 K in applied fields up to 7 T. Heat capacity 
was studied by a relaxation technique (Quantum Design 
PPMS) in the temperature range 2 — 100 K in zero mag- 
netic field. 

The electron spin resonance (ESR) measurements were 
performed with a standard continuous-wave spectrom- 
eter at X-band frequencies (y « 9.5 GHz) by using a 
cylindrical resonator in TE012 mode. The ESR spec- 
tra were analyzed by fitting them with Lorentzian lines. 
From these fits, the linewidth A_B (half-width at half- 
maximum), the resonance field Bres (providing the ESR 
(7- factor), and the intensity (given by the area oc AAB^ 
under the ESR absorption, where A is the amplitude of 
the Lorentzian) were determined. 

Zero- and longitudinal-field (ZF&LF) /iSR experi- 
ments were performed at the ttMB beam line at the Paul 
Scherrer Institute on the GPS spectrometer (Villigen, 
Switzerland). Forward, backward, up, and down positron 



detectors were used for monitoring the asymmetry sig- 
nal A{t). The ZF measurements were performed in the 
transverse mode with muon spin perpendicular to its mo- 
mentum, while the LF measurements were done in the 
LF mode with the muon spin parallel to its momentum. 
The ZF spectra were obtained in the temperature range 
of 5 — 200 K. while the LF measurements were done in 
a series of fields ranging from 1 mT to 640 mT at 100 
K and at 30 mT in the temperature range of 5 to 250 
K. Typical counting statistics were ^ 5 x 10^ positron 
events per each particular data point. 

The ground state of the proposed microscopic mag- 
netic model was studied by coupled-cluster and Lanc- 
zos diagonalization calculations using the program pack- 
ages "The crystallographic CCM" (by D.J.J. Farnell and 
J. Schulenburg) and J. Schulenburg's spinpack, respec- 
tively. To evaluate the Neel temperature, we simpli- 
fied our spin model (see Sec. IIV Ap and performed quan- 
tum Monte Carlo simulations using the loop algorithro22, 
of the ALPS package.— Neutron-scattering patterns for 
possible magnetic structures of CuNCN were calculated 
with FullProf ^ 
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FIG. 2. (Color online) Magnetic susceptibility of CuNCN 
measured up to 650 K in the applied field of 7 T (open circles) 
and the fit with Eq. ((TJ (solid line) along with individual 
contributions described in the text (dashed and dotted lines) . 



III. EXPERIMENTAL RESULTS 
A. Magnetic properties 

Liu et alj^ reported magnetic susceptibility (x) mea- 
surements for CuNCN up to 320 K. They observed 
a nearly temperature-independent susceptibility above 
100 K, whereas at lower temperatures the data showed a 
bend around 70 K followed by a Curie-like paramagnetic 
upturn. Additionally, the increase in the magnetic field 
induced a systematic reduction in Xi thereby indicating 
a minor ferromagnetic contribution of unknown origin. 
We extended the measurements by Liu et al^ in two 
aspects. First, we collected the data above 320 K, since 
there is strong evidence for large couplings on the order of 
2000 Kfi^ so that an evaluation of exchange parameters 
would require the high-temperature data. Second, we 
performed the measurements for several samples to sort 
out external contributions and to better understand the 
role of impurities. We start with the high-temperature 
data, while the sample dependence will be discussed in 
Sec.lnTBl 

Fig. [2] presents the susceptibility data collected up to 
650 K, which is close to the decomposition tempera- 
ture of CuNCNj^ To check for the stability of the sam- 
ple in this temperature range, we performed the mag- 
netization measurements both on heating and on cool- 
ing. Both datascts perfectly matched, thus confirm- 
ing the intrinsic nature of the observed signal. Above 
250 K, the magnetic susceptibility of CuNCN increases. 
This contradicts the itinerant scenario requiring an es- 
sentially temperature-independent Pauli paramagnetism, 
and challenges the conclusion of Ref. [ij on the non- 
magnetic nature of Cu+^ cations in CuNCN. The in- 



crease in X is a clear signature of localized spins that 
are subject to strong AFM couplings inducing a low- 
dimensional and/or frustrated behavior. Quantum fiuc- 
tuations in low-dimensional/frustrated magnets typically 
lead to a broad susceptibility maximum. In our case, the 
maximum lies above the accessible temperature range^^ 
(Fig. [5]), and the measurement can only probe the sus- 
ceptibility below the maximum. Our measurements are, 
therefore, high-temperature on the usual experimental 
scale, but cover a relatively narrow temperature range on 
the scale of the leading exchange coupling J ~ 2500 K 
{T/J < 0.3). 

Although the data in the narrow temperature range 
can hardly be used to establish details of the magnetic 
model, it is still possible to fit the experimental suscepti- 
bility using a given model, and evaluate the relevant ex- 
change couplings. Since DFT calculations^^ consistently 
find the leading AFM coupling J ~ 2500 K along the c di- 
rection, we use the model of a uniform Heisenberg spin-i 
chain supplemented by a temperature-independent con- 
tribution xo for core diamagnetism/van Vleck paramag- 
netism as well as a Curie contribution C /T for the low- 
temperature upturn related to the non-interacting impu- 
rity spins: 



X 



Xo + ^ 



C , A^ASVI 



J 



Xchain(r/J). 



(1) 



Here, g is the g-factor, Na is Avogadro's number, ^b 
is Bohr magneton, and Xchai-a{T/J) is the susceptibil- 
ity of a uniform Heisenberg chain, as given in Ref. l24l . 
Our fit yields xo = 7.3 x 10"^ emu/mol, g = 2.2, 
C = 0.0021 emu K/mol, and J = 2310 K. The fitted g- 
value is in good agreement with the ESR estimate g ~ 2.1 
(see Sec. IIIIDJ) .^^ whereas our experimental estimate of 
J almost perfectly matches the computational prediction 




40 60 

Temperature (K) 

FIG. 3. (Color online) Specific heat of CuNCN measured in 
zero magnetic field. 
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J ~ 2500 K of Ref. |l5|. This result provides a strong 
experimental support for the proposed microscopic mag- 
netic model. 

An important consequence of our model is the low- 
temperature LRO driven by interchain couplings (see also 
Sec. Iivp j^ The bend of the susceptibility curve around 
70 K (Fig. SJ upper panel) is a possible signature of such 
ordering. However, the intrinsic nature of the suscepti- 
bility anomaly has been questioned by Liu et al.^ be- 
cause any visible effect in the specific heat (Cp) is lacking. 
Fig. [3] shows that the heat capacity of CuNCN is indeed 
smooth up to 100 K, but this observation does not ev- 
idence against the LRO. In Sec. HI we mentioned that 
the ordering transitions in quantum magnets may be in- 
visible for thermodynamic measurements because of the 
diminutively small entropy release at the low transition 
temperatures T^ I J — 0.03 (see also Sec. IIVA]) . This is 
the case for CuNCN, where local probes give clear indi- 
cations for the formation of static magnetic fields below 
70 - 80 K (see Sections HHC] and HIlD]). 



B. Sample characterization 

Prior to discussing the /iSR and ESR results, we 
will comment on the problem of possible impurities in 
the polycrystalline samples of CuNCN. The comparative 
study of several samples demonstrated that the magni- 
tude of the low-temperature Curie-like impurity contri- 
bution is strongly sample-dependent. For a further study, 
we selected samples 1 and 2 showing the least pronounced 
and most pronounced low-temperature upturns, respec- 
tively (Fig. 21). Fitting the low-temperature susceptibility 
with the Curie law, we roughly estimate the amount of 
an effective spin-^ impurity as 0.1 % for sample 1 and 
1.3 % for sample 2. It is worth noting that the bend 
around 70 K, which is sometimes thought to be of im- 
purity origin^i^iH is in fact more pronounced in sample 
1 with the significantly lower impurity contribution. In 
sample 2, the anomaly is masked by the stronger impu- 



FIG. 4. (Color online) Magnetic susceptibility for two sam- 
ples of CuNCN with different chemical composition (Table |l|, 
different microstructure (Fig. [S]), and different concentration 
of paramagnetic impurities (see text for details) . The dashed 
line in the upper panel denotes the minimum susceptibility 
and underscores the weak increase in x above 300 K (com- 
pare to the high-temperature data in Fig. [2]). The circle marks 
the bend of the curve around 70 K. 



TABLE I. Chemical composition (in wt.%) of the CuNCN 
samples. 





Cu 


N 


C 


H 
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Sample 1 


60.8 


27.1 


11.6 


0.1 


0.0 


Sample 2 


59.4 


24.5 


11.1 


0.4 


4.7 


Nominal 


61.4 


27.0 


11.6 


0.0 


0.0 



rity signal. 

The powder XRD suggests that both samples 1 and 
2 are single-phase, although this method is sensitive to 
crystalline phases, only. To test the bulk composition of 
the samples, we used the chemical analysis. Represen- 
tative results are reported in Table HI Sample 1 is close 
to the nominal composition, whereas sample 2 consider- 
ably deviates from the anticipated CuNCN composition 
and contains appreciable amounts of oxygen and hydro- 
gen. While the deviations in the composition of sample 2 
could be ascribed to an amorphous impurity phase invis- 
ible for XRD, the actual situation is more complex. The 
two samples not only differ in the chemical composition, 
but also contain CuNCN-type phases with different fea- 
tures of the microstructure. The differences in the crys- 
talline components of samples 1 and 2 are evidenced by 
high-resolution XRD, and further confirmed by the /xSR 
and ESR studies (Sections IIII CI and IIII Dl respectively). 

To better characterize the crystalline CuNCN phase in 
samples 1 and 2, we performed XRD measurements using 
a synchrotron source that provides an excellent resolu- 
tion and shows remarkable sensitivity to the microstruc- 
ture. The room-temperature powder patterns of samples 
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FIG. 5. (Color online) Synchrotron XRD patterns for two 
samples of CuNCN with different chemical composition. In 
the left panel, the pattern of sample 1 is offset for +5000 for 
clarity. Sample 1 shows a pronounced anisotropy of reflec- 
tion halfwidth (note the 004 and 112 reflections in the left 
panel), whereas the reflections of sample 2 demonstrate more 
isotropic angular dependence of the halfwidth and peculiar 
reflection shape (right panel). 



1 and 2 are shown in Fig. [S] In sample 1, we observed 
a dramatic anisotropy of the reflection halfwidth. The 
reflections hkl with both k and I non-zero arc systemat- 
ically broadened, as seen, e.g., from a comparison of the 
two neighboring reflections, 004 and 112 at 29 ~ 9.7°, or 
021 and 110 at lower 26. Reflections of sample 2 look 
more isotropic in terms of the halfwidth, yet their shape 
is somewhat unusual (see the left panel of Fig. [5]) and 
precludes the conventional Rietveld refinement. We find 
that both samples 1 and 2 are not perfectly crystallized 
and exhibit peculiar features of the microstructure. 

At this point, we also comment on the sample charac- 
terization reported in previous studies of CuNCN jiiii^iii 
In sample 1, the chemical composition as well as the mag- 
nitude of the Curie-like susce ptib ility upturn are compat- 
ible with the data of Refs.[lllll2l. andllTl thereby indicat- 
ing similar sample quality. Our sample 2 may be thought 
to be of lower quality, although the deviations in its chem- 
ical composition coexist with - and probably signal - a 
peculiar change in the microstructure of the crystalline 
CuNCN phase. We argue that the previously overlooked 
problem of defects and microstructure might be of crucial 
importance for understanding the magnetism. Although 
the authors of Refs. [ill and [13 describe their samples 
of CuNCN as having "excellent crystallinity" evidenced 
by noticable high-angle Bragg peaks in laboratory x-ray 
data obtained with CuKqi radiation, this test is not sen- 
sitive to extended defects observed in our study. Similar 
to Ref. [ll|, the synchrotron data for samples 1 and 2 
show a number of sharp reflections up to g ~ 7.5 A~^ 
{29 ~ 133° for CuKqi radiation), but this neither pre- 
cludes the dramatic broadening of the low-angle reflec- 
tions with non-zero k and I nor excludes the peculiar 
reflection shape observed in sample 2. 

The variable crystallinity of CuNCN is unrelated to the 
particle size, which would lead to a uniform broadening 
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FIG. 6. (Color online) ZF fiSK time spectra of sample 1 
measured at 5, 65, 90, and 180 K. The solid lines are fits to 
the data with Eq. ^. 

of all Bragg peaks, but rather concerns specific defects re- 
solvable in a high- resolution XRD experiment, only For 
example, the powder XRD data reported in Refs. [ill and 
[TtI show the 004 and 112 reflections at 26' ~ 9.7° as a sin- 
gle peak, and do not allow to observe the difference in the 
halfwidths of these reflections. Therefore, a further char- 
acterization with high-resolution XRD is desirable. We 
believe that genuine "powder samples of CuNCN with 
excellent crystallinity" remain a challenging topic, and 
the problem of microstructure requires further careful 
consideration that lies beyond the scope of the present 
study. 

C. /iSR spectroscopy 

Figure [6| shows ZF /iSR time spectra of sample 1 mea- 
sured at 5, 65, 90, and 180 K. At 5 K, a strongly depo- 
larized /xSR signal is evident, suggesting a magnetically 
ordered state of the sample. With increasing tempera- 
ture, the muon depolarization gradually vanishes, and a 
low relaxation rate at 180 K identifies a complete tran- 
sition to the paramagnetic state. The time dependence 
of the ZF /iSR time spectra is well described with the 
following equation;^ 



A{t) =Ai 



2 1 

- exp(-ATO cos(7^Bint<) + - exp(-ALi) 



y42exp(-A2t). 



(2) 
(3) 



The first term gives the signal of the magnetically- 
ordered polycristalline sample, with 2/3 oscillating and 
1/3 non-oscillating parts, while the second term de- 
scribes the paramagnetic fraction of the sample. Initial 
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FIG. 7. (Color online) Fourier transform (FT) of the oscil- 
lating (2/3) fraction of the /iSR time spectra measured at 5, 
20, 35, 50, and 70 K. The solid lines are FT of the oscillating 
fraction from fits to the data with Eq. ([2]). The non-zero first 
moment of the internal field distribution is evident at 5 K and 
gradually vanishes above 50 K. For better visualization, the 
FT spectra below 70 K are shifted vertically. 



asymmetries Ai and A2 are proportional to the volume 
fractions of the magnetically-ordered and paramagnetic 
phases, respectively, and add up to the temperature- 
independent total asymmetry Ai + A2 — 0.21. The 
transverse relaxation rate At is proportional to the dis- 
persion of the internal magnetic field Bint ■ Al is the lon- 
gitudinal relaxation rate of muon polarization, whereas 
7^ = 27r • 135.53 MHz/T is the muon gyromagnetic ratio. 
The parameter A2 describes the muon relaxation rate in 
the paramagnetic fraction of the sample. It is related to 
the static and dynamic local fields<21 

In contrast to the recent study by Zorko et al.^ we 
clearly identify a nonzero internal field at 5 K in sample 
1. Fig. [7] shows the Fourier transform (FT) for the oscil- 
lating fraction (2/3) of the /iSR-time spectrum (i.e., the 
paramagnetic fraction and the non-oscillating (1/3) frac- 
tion of the signal were subtracted from the full ^SR time- 
spectrum). The non-zero internal field is a direct proof 
for the antiferromagnetically ordered rather than a spin- 
glass state of stoichiometric CuNCN. We note, however, 
the broad distribution of internal fields that may be re- 
lated to several inequivalent muon stop sites in the struc- 
ture and/or to a structural inhomogeneity. This problem 
is further discussed in Sec. IVl 

Figure m shows the ZF /xSR results for Bint, At, Al, 
and Al evaluated using Eq. ([2]). The relaxation rate 
of the paramagnetic fraction A2 is zero below 60 K and 
very close to Al above 60 K. Therefore, we set A2 ~ Al 
above 60 K to avoid correlations between these param- 



eters. At 5 K, 100 % volume fraction of the sample is 
in the AF ordered state (see Ai{T) in the left panel of 
Fig. |5|). With increasing the temperature above 25 K, 
the magnetic volume fraction gradually decreases to 50% 
at about 62 K and to ~ 25% at 70 K. At higher tem- 
peratures, Al decreases linearly and eventually vanishes 
at 180 K. By contrast, the internal magnetic field i?int 
probed at the muon site(s) decreases with temperature 
from 30 mT at 5 K to at ~ 70 K. Above 70 K, spins 
are in a highly disordered state, as seen from the tem- 
perature dependence of the transverse relaxation rate At 
(right panel of Fig. |S]). The longitudinal relaxation rate 
also shows a small peak at about 50 K, as expected at the 
proximity to a magnetic transition. This peak is smeared 
due to a rather broad magnetic transition. Presently, we 
are unable to find a unique reason for the marginal frac- 
tion of the ordered phase (non-zero ^1) between 70 and 
180 K. Plausible explanations would be the formation of 
a muon-induced magnetic phase in small regions of the 
sample or the structure inhomogeneity. The zero inter- 
nal field as well as the well-defined ESR absorption above 
70 K (Sec, nil D|) corroborate the lack of the macroscopic 
LRO above 70 K. 

The /iSR time-spectrum of the paramagnetic fraction 
is described with an exponential function [see Eq. ^]. 
The exponential relaxation rate is characteristic either 
for the dynamic muon relaxation or for the relaxation due 
to static diluted magnetic centers. In order to distinguish 
between these two possibilities, we performed the LF /iSR 
experiment (also referred as a decoupling experiment) at 
100 K where the dominant fraction of the sample is in the 
paramagnetic state. The /iSR time spectra were analyzed 
with the equation: 



yl(t) = yl'exp(-Ai). 



(4) 



The fit was performed for the times above 0.5 /zs to ex- 
clude the influence of the fast relaxing component of the 
/iSR time-spectrum (see Fig. [5]). The longitudinal-field 
(Bl) dependence of A is shown in the inset of Fig. [HI 
The /iSR signal is decoupled at about 0.5 mT suggesting 
the dominant static character of muon depolarization in 
the paramagnetic state. Therefore, to study the weak 
dynamic muon relaxation we performed a temperature 
scan in the applied longitudinal field B\^ = 30 mT where 
static internal fields are completely decoupled. The fit- 
ting results are summarized in Fig. [S] The maximum 
relaxation rate A is 0.02 /is~^ which is close to the in- 
strument resolution limit resulting in substantial error. 
The temperature dependence of A has a peak at about 
50 K, similar to one observed on Al in ZF /iSR (right 
panel of Fig. [8]). 

We also measured the ZF /iSR spectrum for the off- 
stoichiometric sample 2 at 5 K (Fig. [TU)) . In contrast to 
sample 1 , sample 2 retains a weak paramagnetic compo- 
nent down to low temperatures, and shows a more shal- 
low dip in A{t). The latter feature indicates a broader 
distribution of internal fields, as confirmed by the Fourier 
transforms shown in the insets of Fig. fTUl In sample 1, the 
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FIG. 8. (Color online) Summary of ZF fiSK results fitted with Eq. 0. Left panel: temperature dependence of magnetic volume 
fraction Ai and internal magnetic field Bint- Right panel: longitudinal (Al) and transverse (At) relaxation rates as a function 
of temperature. 
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rate A as obtained with Eq. (g} in the field of Bl = 30 mT. 
The inset shows the decoupling LF experiment at 100 K. The 
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FIG. 10. (Color online) ZF /iSR spectra of samples 1 and 2 
measured at 5 K and the respective Fourier transforms show- 
ing internal fields -Bint. 



D. ESR 



distribution of Bi^t has a well-resolved maximuni around 
30 mT, whereas in sample 2 internal fields between and 
20 niT are observed with nearly equal probabilities. The 
different /iSR spectra confirm that the primary differ- 
ences between the samples are unrelated to the possible 
contamination with an amorphous paramagnetic impu- 
rity, but rather reflect the differences in the crystalline 
"CuNCN" phase. 



While muon scattering directly probes the amount of 
the magnetically-ordered phase, this experiment involves 
the interaction between the crystalline material and im- 
planted muons. Therefore, extrinsic phenomena driven 
by muons and, particularly, the local magnetic order in- 
duced by muons can not be excluded. To check for the 
intrinsic nature of the observed magnetic transition, we 
performed the ESR measurements. 

The ESR spectra of the two samples of CuNCN 
(Fig.llip show clear differences except for the narrow line 
with the temperature-independent linewidth of 20 niT at 
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FIG. 11. (Color online) X-band ESR results of CuNCN (sample 1, left panels and sample 2, right panels). The ESR signals 
(field derivative of absorbed microwave power) were fitted with a Lorentzian shape (dashed lines) using a two lines (sample 1) 
and three lines (sample 2), as further explained in the text. The lower panels display the corresponding parameters linewidth 
(AS), g-factor (determined by the resonance field) and intensity (integrated microwave absorption). 



a field of 327 mT (5 = 2.04(2)). These parameters are 
typically expected for the resonance of non-interacting 
Cu+2 spins. The left panel of Fig. HU shows the ESR 
spectrum of sample 1 which consists, beside the narrow 
line, of one broad Lorentzian line denoted by the dashed 
curve in the experimental spectrum. The parameters of 
this broad line display a temperature dependence char- 
acteristic to the presence of the low-temperature mag- 
netic order. The pronounced increase (decrease) in the 
linewidth ((^-factor) indicates the slowing down of spin 
fluctuations and the onset of internal magnetic fields in 
agreement with the //SR observations. Upon cooling, the 
line intensity is slowly decreasing and eventually drops 
down to zero below 70 K. Therefore, the intrinsic mag- 
netic susceptibility of CuNCN vanishes at low temper- 
atures, as typical for the LRO state. Another possible 
interpretation is the presence of a spin gap that leads to 
the decrease in the magnetic susceptibility at low tem- 
peratures and, consequently, to the vanishing ESR line. 
Overall, the temperature evolution of the two lines in the 
ESR spectrum is consistent with the magnetic suscepti- 



bility of sample 1 (Fig. HI upper panel). Based on the 
observation of two different lines, we are able to separate 
the experimental susceptibility into the Curie-like para- 
magnetic contribution and the intrinsic contribution of 
CuNCN that vanishes below 70 K. 

The ESR spectra of sample 2 could be fitted by three 
Lorentzians as shown by the solid lines in the upper right 
panel of Fig. [TTJ Line 3 is the aforementioned narrow 
temperature-independent line. This line is superimposed 
by a strong and broad line 2 dominating the spectra 
in the whole investigated temperature range. In agree- 
ment with the larger susceptibility of sample 2 (Fig. U]) , 
this line causes an order of magnitude larger ESR inten- 
sity for sample 2 compared to sample 1. Despite similar 
linewidths, line 2 in sample 2 is rather different from the 
broad line in sample 1. While the latter disappears below 
70 K, the former is observed down to 5 K, with marginal 
changes in the linewidth and g-value. Interestingly, the 
ESR spectrum of sample 2 features an additional contri- 
bution (line 1) that does vanish below 70 K, yet showing 
a smaller linewidth and smaller resonance field (higher g- 



factor) than the intrinsic signal (broad hnc) in sample 1. 
Therefore, the temperature of 70 K is also characteristic 
for sample 2, although the magnetic transition (or the 
spin gap) takes place in a part of the sample only. This 
compares well to the incomplete and rather inhomoge- 
neous magnetic ordering observed with /iSR (Fig.lTOland 
Sec lmCl) . 

Presently, we are unable to decide whether lines 1 and 
2 indicate different magnetic phases of sample 2, or ev- 
idence some peculiar magnetism within the single off- 
stoichiometric CuNCN phase. However, both ESR and 
/iSR prove that sample 2 can not be simply considered as 
the crystalline CuNCN phase mixed with an amorphous 
impurity. The deviations in the chemical composition 
and/or microstructure have strong influence on the mag- 
netism of CuNCN. 



IV. THEORETICAL ANALYSIS 

The ESR and /iSR results give strong evidence for 
the emergence of static magnetic fields in CuNCN below 
70 K. This onset temperature matches the weak suscepti- 
bility anomaly. However, the specific heat and neutron- 
scattering dat a^^'^'^ do not show any signatures of the 
LRO down to 2 K. In the following, we will apply a va- 
riety of numerical techniques to investigate the ground 
state and thermodynamic properties of CuNCN. Subse- 
quently, we use the results of this analysis to address the 
sensitivity of different methods to the anticipated long- 
range magnetic order, and suggest a consistent interpre- 
tation for the available experimental data. 



A. Ground state 

According to Ref. |T^ the microscopic magnetic model 
of CuNCN entails AFM spin chains along the c direc- 
tion. The chains arc coupled by the nearest-neighbor 
ferromagnetic (FM) interaction Ji ~ —500 K and the 
next-nearest-neighbor AFM interaction J2 ~ 100 K along 
a. Additionally, the diagonal coupling Jac — J2 con- 
nects the chains in the ac plane, whereas a very weak 
AFM coupling Jf, ~ 5 K links the neighboring layers in 
a frustrated manner (see the bottom panels of Fig. [T]). 
The spin lattice is quasi-one-dimensional, yet showing a 
strong anisotropy of the interchain couplings along the a 
and h directions. Since Jh is at least an order of magni- 
tude smaller than the couplings along a, we start with 
analyzing the two-dimensional (2D) spin lattice in the ac 
plane. 

The couplings J, Ji , and Jac favor the columnar AFM 
order featuring parallel spins along a and antiparallel 
spins along c (Fig. [Ij. However, the columnar order is 
frustrated by J2 that prefers antiparallel arrangement of 
next-nearest-neighbor spins along a. To find the classical 
ground state of this magnetic model, we write down the 
energy as a function of fc^ and kz , the components of an 



arbitrary propagation vector k = {k^, kz) in the ac plane. 
The energy per the unit cell of the spin lattice is: 



E = -[J cos kz + Ji coskx + J2Cos(2fc3;)-|- 

+ 2JacC0Sfc2:C0Sfcz)]. 



(5) 



The energy minimum is at k = (0, tt) for 
•h/i-Ji + 2Jqc) < J (columnar AFM phase) 
and at k = [arccos [( — Ji + 2Jac)/4J2] , tt] for 
J2/{—Ji + 2Jqc) > \ (spiral phase). Comparing 

this result to the better studied model of coupled 
frustrated spin chains (FM Ji, AFM J2, AFM J, and 
Jac = 0))^ we note a similar competition between 
columnar (coUinear) and spiral phases that are sepa- 
rated by a quantum critical point at J2/J1 = — 4. The 
diagonal coupling Jac shifts the quantum critical point 
to Jil{.—J\ + 2 Jac) = \ without changing the nature of 
the competing phases. Note that the frustration is only 
present along the a direction, hence the leading coupling 
J along c does not influence the stability of the collinear 
and spiral phases. 

The spin lattice of CuNCN with J2/(- Ji + 2Jac) "^ 
0.14 clearly favors the columnar AFM phase. While the 
classical model suggests the ordered moment of 1 /i^ in- 
dependent of individual exchange couplings, quantum ef- 
fects have to be taken into account for a spin-i system. 
To ensure the reliable treatment of such effects, we use 
the coupled-cluster method (CCM). An efficient applica- 
tion of CCM to low-dimensional spin systems has been 
illustrated in Refs. I28l436l . Here we give only an outline 
of the main relevant features of the CCM. 

The core of the CCM method is the choice of a nor- 
malized reference state |$) together with a complete set 
of (mutually commuting) multi-configurational creation 
operators {C^} and the corresponding set of their Her- 
mitian adjoints {Cl\- In our case, the reference state |$) 
is the columnar AFM phase of the classical spin model 
(parallel spins along a, antiparallel spins along c). Then 
we perform a rotation of the local axis of the spins such 
that all spins in the reference state are directed along 
the negative z axis, i.e., in the rotated coordinate frame 
the reference state reads |$) = | J,) | J,) | J,) The corre- 
sponding multispin creation operators then can be writ- 
ten as C/ = s+, s^st^ s'^ststf ,■ ■ ■ ^ where the indices 
a, /?, 7 . . . label arbitrary lattice sites. 

The kct- and bra- ground states are given by: 



|vl') = e^|$), S = Y,SjC] 
(§1 = ($|5e-^, 



/#o 



S^\ + ^SiCi. (6) 



Using ($|C/ = = C/l*) V/ 7^ 0, Cp^ = 1, the com- 
mutation rules [C^,C^] = = \Cl,Ck\ the orthonor- 



mality condition ($|C/Cj |<I> 



and completeness 



/ 



C+|$)($|C7 = 1 = |$)($| +^C+|$)($|C/, we get 
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FIG. 12. (Color online) COM LSUBn results as well as ex- 
trapolated (n -^ oo) data for the ordered moment m at 
Ji/J — —0.2. The couplings J2 = Jac are given in units 
of J = 1. The two extrapolation schemes are explained in the 
text. The shaded bar denotes the likely range for the ordered 
magnetic moment in CuNCN. 



a set of non-linear and linear equations for the eorrelation 
coefficients Sj and Sj, respectively. The order parame- 
ter (sublattice magnetization or ordered moment, m) is 
given by 



m 



1 ^ 
N ^^ 



:*i<j*> 



(7) 



where sf^ is the spin operator expressed in the rotated 
coordinate system, and (i, n) denotes the position of the 
lattice site. The CCM provides results for infinite lat- 
tices. The only approximation of the CCM is the trun- 
cation of the expansion for the correlation operators S 
and S. We use the well-established LSUBn scheme that 
includes all multispin correlations on the lattice with n 
or fewer neighboring sites. To account for all multispin 
correlations, the CCM estimates obtained on different 
levels of the LSUBn approximations (to„) are then ex- 
trapolated to the n — )■ 00 limit (rrioo) using the extrap- 
olation schemes m,i = rrioo + a/n + b/n^ (scheme 1) 
and nin ~ nioo + a/n^^^ + b/n?'^'^ (scheme 2). The for- 
mer scheme is typically used for quantum spin systems 
with the AFM long-range oidei^— whereas the latter 
scheme has been successfully applied to systems near a 
quantum critical point, where the magnetic order param- 
eter m is small^i^i^ 

Fig. [12] shows the extrapolated ordered moment m — 
rUoo for Ji/J ~ —0.2 and variable J2 = Jac- We find 
that the frustrating coupling J2 has no appreciable ef- 
fect on the ordered moment m in the relevant parameter 
range. Depending on the extrapolation scheme used, we 
find m = 0.34 — 0.46 /is, where the largest value is ob- 
tained for scheme 1 at J2 = Jac = and the lowest value 
is found for scheme 2 at J2 = Jac = 0.5Ji. These values 
are well below the classical value of 1 ^b and also signifi- 
cantly lower than the ordered moment of ■~0.6 fiB on the 
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FIG. 13. (Color online) CCM-LSUB8 and ED resuhs for the 
spin-spin correlations (sqSr) along the c (top panel) and a 
(bottom panel) directions at Ji/J — —0.2. The couplings 
J2 — Jac are given in units of J = 1. The ED results for 
R = (0, 4) are not available because of the limited size of the 
finite lattice. 



spin-i square lattice (see, for example, Refs. [30, [S^, and 



spi i 



371 ). This small value of the order parameter can be at- 
tributed to the dominating exchange coupling J along the 
c direction that leads to a quasi-one-dimensional system 
of weakly coupled antiferromagnetic spin- 2 Hcisenberg 
chains. Simple models of such weakly coupled chains are 
well studied, e.g., in Refs. l34l I38l - l4ll . These works show 
that even an infinitesimally small (non-frustrated) inter- 
chain coupling leads to the magnetic long-range order at 
zero temperature. Although our model does not corre- 
spond directly to these simple models considered in the 
previous literature, a similar conclusion is also valid for 
our model, in particular, since our numerical simulations 
indicate that the effect of the frustrating coupling J2 is 
very weak. 

To get additional insight into the magnetic ordering 
and the role of the frustrated coupling J2, we present 
spin-spin correlation functions (sqSr) for various lattice 
vectors R (Fig. [T^ . In addition to the CCM data, we 
also show the results of Lanczos exact diagonalization 
(ED) for a 8 X 4 finite lattice with periodic boundary 
conditions.— Obviously, the spin-spin correlations sup- 
port the anticipated columnar AFM ground state with 
a large negative correlation for nearest neighbors along 
c and a sizable positive correlation for nearest neighbors 
along a. Most correlations are insensitive to J2 and Jac, 
although the correlation on the Ji bond is slightly en- 
hanced as J2 and Jac increase. This signifies a slight en- 
hancement in the coupling along the a direction, yet this 
effect is countered by the increasing frustration. Over- 
all, the weak changes in the spin-spin correlations are 
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consistent with the nearly constant value of the ordered 
moment, as evaluated by CCM (Fig. [T2). Note that the 
weak correlations along the a direction support the pic- 
ture of weakly-coupled spin chains leading to a small or- 
dered moment. Although the ED and CCM results for 
spin-spin correlations show very similar qualitative be- 
havior, the quantitative differences are obvious. These 
differences can be attributed to strong finite-size effects 
in the ED data, since in the c direction we have a spin 
ring of 4 sites, only. 

Using the CCM and exact diagonalization results, we 
conclude that the spin lattice of CuNCN develops the 
columnar AFM order in the ac plane. This order cor- 
responds to the propagation vector k ~ (0,fcj,,7r) with 
respect to the unit cell of the spin lattice (Fig. [II bottom 
right panel) or k = (0, ky, 0) with respect to the unit cell 
of the atomic structure. Now, we comment on the possi- 
ble interlayer ordering along the b direction. According 
to the left panels of Fig. [I] the weak interlayer coupling Ji, 
forms isosceles triangles and therefore frustrates the in- 
terlayer order on the classical level. This classical degen- 
eracy is lifted by quantum fluctuations favoring coUinear 
ground states via the order-from-disorder mechanismi^ 
We thus expect the collinear order in CuNCN, although 
the periodicity along b cannot be determined unambigu- 
ously. Note that a similar case of weak and frustrated 
interlayer couplings has been reported for the Heiscnberg 
model on the body-centered tetragonal lattice.— An elab- 
orate spin-wave study of Ref. |4J suggests two collinear 
states with the doubled and quadrupled periodicity along 
the interlayer direction. Such states are strongly favored 
over any non-collinear configurations, yet the preference 
for the state with the quadrupled periodicity appears in 
high orders of spin-wave theory, only. Therefore, we ex- 
pect the doubled or quadrupled periodicity along the b 
direction in CuNCN, i.e., propagation vectors are k = 
or k = (0, TT, 0) with respect to the crystallographic unit 
cell of CuNCN. 



B. Thermodynamic properties 

CCM is a handy tool for studying frustrated spin 
systems, but its applications are restricted to ground- 
state properties. To investigate the finite-temperature 
behavior of our microscopic model, we have to apply 
ED that puts severe restrictions on the accessible lattice 
size (note, for example, the pronounced finite-size effects 
for the spin-spin correlations in Fig. I13p . However, the 
ground-state properties are weakly dependent on J2 and 
Jac- Therefore, we can safely simplify our model by re- 
moving the coupling J2 and eliminating the frustration. 
This simplification enables us to apply quantum Monte- 
Carlo (QMC) techniques that readily treat large systems, 
thus effectively overcoming finite-size effects. Since the 
2D model in the ac plane maintains the long-range order 
at zero temperature only, we also introduce a realistic 
interlayer coupling J^^/J — 0.002 that should stabilize 
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FIG. 14. (Color online) Binder cumulant (rn^) / {rrvi)'^ calcu- 
lated for different system sizes at Ji/J = —0.2, Jac/ J ~ 0.04, 
and the effective interlayer coupling J^ /J — 0.002. The 
shaded bar shows the estimated Tn/J — 0.12. 



the LRO above K. 

The spin lattice of CuNCN is spatially anisotropic, 
with the leading coupling along c, the five times weaker 
coupling along a, and a diminutively small coupling 
along b. Following Refs. |40, |4a, and |4a, we adjust 
the aspect ratio of finite lattices to account for differ- 
ent correlation lengths along a, 6, and c. Specifically, 
we used L/4 x L/16 x L lattices with L < 160 (up to 
64000 sites). The Neel temperature T/v is determined 
from the standard scaling procedure for the Binder ratio 
B = (m^)/(?Ti^)^, where iris is the staggered magneti- 
zation. Since B{T) is independent of the lattice size at 
T = T/v,— the Neel temperature is precisely determined 
as the crossing point of the B{T) curves calculated for 
different L. Fig. [T3] shows the scaling procedure and 
the resulting Tfq / J ~ 0.12. Using the experimental es- 
timate of J ~ 2300 K, we arrive at Tn — 275 K, which 
is much larger than the experimental value of 70 K, be- 
cause the non-frustrated interlayer coupling Jf, is consid- 
ered. In CuNCN, the frustrated couplings Jh will impede 
the LRO, thus reducing T/v and bringing it closer to the 
experimental value. 

It is instructive to compare the above numerical result 
to previous theoretical studies of coupled spin chains 3 
Despite the abundance of spatially anisotropic interchain 
couplings in model quasi-lD compounds;^^— theoreti- 
cal works rather focus on the case of isotropic interchain 
couplings (same coupling along a and &) ,40,52,53 because 
systems with large spatial anisotropy are quite difficult 
to model. To facilitate the comparison, we set the di- 
agonal coupling Jac to zero and consider a simplistic 
model of uniform chains coupled with Ji/ J = —0.2 and 
J^^/J = 0.002 along a and b, respectively. The resulting 
Tn/J is 0.09 compared to 0.12 in the model including Jac, 
and shows the non-negligible role of Jac in enhancing the 
interchain coupling within the ac plane. We now com- 
pare this estimate to the results of Ref. \5^ for isotropic 
interchain couplings J±. Using the averaged coupling 
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FIG. 15. (Color online) Magnetic specific heat of spin chains 
with a uniform interchain coupling J^/ J of 0.05, 0.1, 0.15, 
and 0.2. The decrease in J^/ J shifts the transition anomaly 
to lower temperatures and dramatically reduces its magni- 
tude. 
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FIG. 16. (Color online) Magnetic specific heat for the spin 
lattice of CuNCN, with Ji/J = -0.2, Jac/J = 0.04, and 
J^ / J — 0.002. The inset shows the magnetic entropy Smag 
obtained by integrating Cms^^/T. 
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tually becomes invisible at sufficiently weak interlayer 
couplings. 

In CuNCN. the very weak coupling along the 6 di- 
rection enhances quantum fluctuations and reduces T/y 
compared to the scenario of the averaged interchain cou- 
pling Jf ^ {\Ji\ + Jb)/'2 ^ O.IJ. This leads to a further 
decrease in the transition anomaly. Fig. [1^] shows the 
simulated specific heat for the relevant exchange param- 
eters of CuNCN, as derived from the DFT calculations j^ 
We were unable to discern any signatures of the tran- 
sition anomaly at T/v/J ~ 0.12. because the magnetic 
entropy available at the Neel temperature is exceedingly 
low. We illustrate this effect in the inset of Fig. [161 where 

T 

/ %a dT is shown. 



0.13 J mol-i K-i is less than 3 % of 
the total magnetic entropy i?ln2 ~ 5.75 J mol~^ K^^. 
Considering the inevitable effect of exchange anisotropy 
(broadening of the transition anomaly) and the large lat- 
tice contribution of about 20 J mol^^ K^^, which is 150 
times larger than S'mag(TV) (Fig. E]), we conclude that 
the specific heat anomaly at the magnetic transition in 
CuNCN can not be observed experimentally, especially 
in experiments on polycrystalline samples. 

Altogether, our QMC simulations are restricted to the 
simplified, non-frustrated spin lattice and may not be suf- 
ficient for the quantitative description of CuNCN. Nev- 
ertheless, we are able to demonstrate that the expected 
Neel temperature is quite low, and the magnetic spe- 
cific heat does not show any signatures of the magnetic 
transition. In CuNCN, the frustrated arrangement of 
the interlayer couplings Jt, will further impede the mag- 
netic order, thus reducing Tjsi and leaving no room for 
the experimental observation of the magnetic transition 
in heat-capacity data. 



C. Neutron scattering 



of J'^ = (|Ji| + Jb)/2 ~ O.IJ, one finds a much larger 
Tjq/ J ~ 0.169.^ This discrepancy is often bypassed by 
assuming an empirical three-fold or four-fold reduction 
in T/v due to the spatial anisotropy of the interchain 
couplings 1^^^ Our numerical results suggest that this 
empirical recipe underestimates T/v /J ~ 0.09, which is 
only twice reduced compared to the estimate from the 
averaged interchain coupling J'^. 

To explore the transition anomaly at Tjq, we evalu- 
ated the magnetic specific heat for coupled spin chains. 
We first consider the isotropic interchain coupling J±^ 
(Fig. [T5|).— At J±/ J = 0.2, the pronounced transition 
anomaly is superimposed on the broad specific heat maxi- 
mum. At lower J^/ J, the T/y is reduced, thereby shifting 
the anomaly below the maximum and shrinking its size. 
This trend follows recent experimental- and theoretical^ 
findings for quasi-two-dimensional spin systems, where 
the transition anomaly is dramatically reduced and even- 



The intensity of the neutron scattering from a mag- 
netic structure is proportional to the square of the mag- 
netic structure factor, which is, in turn, proportional to 
the ordered magnetic moment. Therefore, for low mag- 
netic moments the intensity of magnetic reflections in 
a neutron diffraction experiment decreases dramatically. 
Previous studiesi^ii^ put forward the lack of the mag- 
netic scattering as one of the arguments against the LRO 
in CuNCN. To flnd out whether or not the anticipated 
LRO could be observed in these experimentS )^^'^'^ we sim- 
ulated neutron diffraction patterns of CuNCN. The input 
parameters are the propagation vector and the ordering 
pattern, the ordered moment and its direction, as well 
as the scale factor for the atomic structure. The latter is 
necessary to scale the magnetic reflections with respect to 
the nuclear peaks. The ordering pattern in the ac plane 
and the ordered moment are provided by our theoretical 
study in Sec. IIV Al Since the direction of the magnetic 
moment is presently unknown, different possibilities were 
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FIG. 17. (Color online) Simulated magnetic scattering from CuNCN for the possible propagation vectors k = (left panel) 
and k = (0, i,0) (right panel) and different directions of the magnetic moment. The shaded area denotes the q range probed 
in the polarized-scattering experiment of Ref. [la . The intensity of the largest nuclear reflection is scaled to 1 arb. unit. Note 
that the magnetic scattering is at least 100 times weaker than the nuclear scattering. 



explored. We used the propagation vectors k = and 
k = (0, TT, 0), as explained in Sec. IIVAI 

Simulated patterns for two possible propagation vec- 
tors and three representative directions of the magnetic 
moment are shown in Fig. 1171 The maximum intensity 
of the magnetic scattering is about 1.2 % (hereinafter, 
we measure intensities as fractions of the largest nuclear 
peak) for the magnetic reflection at q = 0.6 — 0.8 A~^ 
(001 + k). At higher q values, the magnetic intensities 
are notably lower because of the rapidly decreasing mag- 
netic form factor of Cu"*"^. However, even the strong 
magnetic reflection at low q may be hard to observe. For 
example, Fig. 1 of Ref. [I^ presents the low-temperature 
neutron-diffraction data for CuNCN. In these data, the 
low-angle 002 nuclear reflection with a sizable intensity 
of 2.3 % is fully masked by the background. This im- 
plies that the data arc by far insufficient to observe the 
magnetic scattering in CuNCN for any of the ordering 
patterns. 

The polarized-scattering data of Ref. [Tj have an addi- 
tional advantage of detecting the magnetic reflections in 
a separate channel, thereby probing the possible mag- 
netic contributions to nuclear peaks and reducing the 
background. An inevitable drawback is the reduction 
in the total neutron flux because of the polarization fil- 
ter. Therefore, the observation of the magnetic scattering 
from a spin-i magnet with the low magnetic moment of 
about 0.4 ^B is at best challenging, at least with the 
instrument used in Ref. [l3l^ We also note that the res- 
olution in the nuclear channel was high enough to detect 
an imperfection in the sample, as evidenced by a for- 



eign peak at g = 0.57 A~^ that brings back the problem 
of structural and chemical homogeneity of large powder 
samples used in neutron experiments. 

In the polarized neutron-scattering experiment, no 
magnetic peaks were observed up to gmax — 2.1 A~-'^j^ 
Even assuming that the experiment of Ref. |l3| were suffi- 
ciently sensitive to detect any magnetic scattering below 
Qmax, we find at least one ordering pattern (k = 0, m 
along c) yielding magnetic refiections above gmax only. 
In this pattern, the lack of the intense reflection at 
q ~ 0.65 A (001 nuclear peak) is because of symmetry re- 
strictions, since the magnetic moment along c forbids the 
magnetic contribution to OOL Therefore, the polarized- 
scattering experiment reported in Ref. Il3| can not be con- 
sidered as an ultimate test for the presence or absence 
of LRO in CuNCN. Neutron diffraction measurements 
for CuNCN largely remain an open problem, and require 
more sensitive instruments accessing larger angle range 
or, preferably, the use of single crystals. 



V. DISCUSSION 

We have shown that CuNCN is a low-dimensional mag- 
net featuring strong quantum fluctuations and an intri- 
cate magnetic transition around 70 K. In the following, 
we propose a plausible microscopic scenario for the mag- 
netism of this compound, and start with analyzing the 
magnetic state of Cu"*"^ . Previous computational studies 
controversially reported a magnetioi^ as well as a non- 
magnetioi^ state for copper. This issue is now resolved 
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by the ESR experiment showing a characteristic absorp- 
tion hne that signals the presence of the locaUzed mag- 
netic moment (Fig. lTl]) . The locahzed behavior is further 
underpinned by the increase in the magnetic suscepti- 
bihty above room temperature (Fig. [2]). Our fit to the 
data suggests a very strong exchange couphng of about 
2300 K. These observations are in remarkable agreement 
with the computational predictions of Ref. llSl and iden- 
tify CuNCN as a charge-transfer insulator that comple- 
ments the broad family of Cu"'"'^-based insulating quan- 
tum magnets. 

A puzzling feature of CuNCN is the magnetic transi- 
tion around 70 K. This transition was first evidenced by 
the bend in the temperature dependence of the magnetic 
susceptibility, but the lack of the corresponding specific 
heat anomaly as well as the missing magnetic scatter- 
ing in the neutron diffraction data led the authors of 
Refs. |lj and |T^ to conclude on the extrinsic nature of 
the susceptibility feature. In Sec. IIV| we have shown 
that neither specific heat nor neutron diffraction mea- 
surements are sufficiently sensitive to the possible LRO 
in CuNCN. Further, our fiSK data reveal the formation 
of static magnetic field in the bulk of the stoichiomet- 
ric CuNCN sample below 70 K. The emergence of static 
fields is accompanied by a vanishing absorption line in the 
ESR spectrum, thus confirming the onset of the LRO. 

Experimental techniques have different sensitivity to 
the LRO magnetic state and may additionally in- 
volve perturbations, such as interactions with implanted 
muons. Therefore, it is important to examine all possible 
scenarios that could explain the experimental data. The 
vanishing ESR absorption is typical for both LRO and 
spin-gap states of a material. While the gapped singlet 
state does not feature any static fields, it is fragile toward 
interactions with muons that show sizable depolarization 
characteristic of a local spin-glass state in the vicinity of 
the implanted muon. Such states arc characterized by 
the zero internal field and the monotonous decrease in 
the muon polarization asymmetry A{t)^ Our data show 
the dip in A{t) that signifies the non-zero internal field 
(Fig. [7]) and oscillations typical for the LRO state. More- 
over, the dynamic relaxation rate A shows the anomaly 
at T^, which is also consistent with the LRO scenario. 
The observation of the non-zero internal field excludes 
the formation of the muon-driven spin glass in CuNCN. 
The spin gap in CuNCN would further contradict the 
microscopic analysis (Sec. IIVI and Ref. Il5l ). because the 
proposed spin lattice leaves no room for the spin gap for- 
mation. 

An opposite situation is known for gaplcss quantum 
magnets, where muons disrupt exchange pathways and 
induce a local magnetic order that results in well-defined 
oscillations of A{t) and few sharp peaks in the Fourier- 
transformed spectrumi^ This is again different from our 
observations (Figs. |6] and [7]). More importantly, the 
muon-induced magnetic order has its own temperature 
scale that is usually different from the intrinsic tempera- 
ture scale of the systemi^ In CuNCN, the onset of static 



magnetic fields, as probed by muons, is accompanied by 
the vanishing ESR absorption. Therefore, the effect is 
clearly intrinsic. At this point, we mention the recent 
results by Zorko et al.^ who also observed the strong 
muon depolarization and the vanishing ESR line below 
70 K. Zorko et alJ^ ambiguously attribute the /xSR sig- 
nal to a spin-glass state induced by the implanted muon, 
although their spectra show a similar dip in A{t) and 
likely lead to the non-zero internal field. Unfortunately, 
the authors of Ref. [ij do not consider the possibility of 
the LRO state, which is the most likely and microscopi- 
cally justified scenario for the ground state of CuNCN. 

We have shown that the available experimental data 
on CuNCN are consistently explained in the framework 
of the microscopic scenario proposed in Ref. llSl . The 
quasi- ID spin lattice leads to strong quantum fluctua- 
tions that reduce the Neel temperature and the ordered 
moment, thus impeding the experimental observation of 
the columnar AFM LRO by heat-capacity measurements 
and neutron scattering. The results of the magnetic sus- 
ceptibility, fiSK, and ESR are a compelling evidence for 
the emergence of static magnetic fields below 70 K, in 
agreement with the developed microscopic scenario. The 
onset of static magnetic fields is further corroborated by 
the broadening of the ^^N nuclear magnetic resonance 
(NMR) line reported by Zorko et al^ 

Despite this robust qualitative scenario, several aspects 
of the low-tcmpcrature magnetic behavior require further 
attention. One obvious problem is the broad distribu- 
tion of internal fields probed by the implanted muons 
(Fig. [7]). This broad distribution is rather unusual for 
the LRO state that typically leads to one or few charac- 
teristic fields at muon stop sites, well-defined oscillations, 
and narrow peaks in the Fourier-transformed spectra. A 
plausible explanation is a large number of inequivalent 
stop sites with different internal fields. For example, the 
/iSR study of spin-i frustrated square lattices in V"*"* 
phosphates^ revealed a similarly broad distribution of 
internal fields in Pb2VO(P04)2 (the variance of the dis- 
tribution is about 10 mT, as compared to 20 — 25 mT in 
CuNCN), although the LRO in this compound is inde- 
pendently confirmed by polarized neutron scattering)^ 
The crystal structure of CuNCN is quite simple, but its 
relation to the possible muon stop sites is far from obvi- 
ous, because the negative charge is associated with rather 
delocalized molecular orbitals of the NCN unit^i^ in con- 
trast to the localized oxygen orbitals in oxides. A fur- 
ther microscopic analysis of the possible muon stop sites 
is highly desirable to achieve more complete understand- 
ing of the /iSR response in CuNCN and other transition- 
metal carbodiimidcs. 

A closely related problem is the broad distribution of 
NMR spin-lattice relaxation times reported in Ref. [l3- 
This finding is based on fitting the magnetization recov- 
ery with a stretched exponent, and taken as an evidence 
for an inhomogeneous magnetic ground state of CuNCN. 
Fitting the magnetization recovery with a stretched ex- 
ponent produces a cumulative parameter j3 that should 
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be equal to 1 for the single relaxation time. The devia- 
tion from unity either describes a broad distribution of 
relaxation times or the presence of several distinct relax- 
ation times in the system under investigation;^ For ex- 
ample, the strongly reduced (3 ~ 0.6 (compare to /3 ~ 0.5 
in CuNCNii) has been observed in the LRO state of 
CaV2 04i^ However, the magnetic ordering in CaV2 04 
is in no way inhomogeneous, as shown by the comprehen- 
sive neutron-scattering studyi^ Therefore, the complex 
behavior of the magnetization recovery in CuNCN could 
also be understood as the presence of several distinct re- 
laxation times in the LRO state. 

The last, and probably most severe problem is the 
smeared magnetic transition in CuNCN. We identify 
T/v ~ 70 K as the onset temperature, although the com- 
plete magnetic order in the bulk of the CuNCN sample 
is established below 20 K, only (Fig.[8l see also Fig. 2 in 
Ref. \Vn \. The broad intermediate region between 20 and 
70 K is somewhat unusual for crystalline magnetic sys- 
tems showing abrupt (and typically second-order) transi- 
tions to the LRO state. The broadening of the transition 
might be related to the peculiar microstructural effects 
observed in our synchrotron XRD experiment (Fig. [5]). 

Altogether, we argue that the bulk of the experimen- 
tal observations on CuNCN can be reconciled within the 
quasi- ID microscopic scenario proposed in Ref. Il5l The 
Cu+^ sites hold localized magnetic moments, and the 
temperature of 70 K is the characteristic temperature 
related to the onset of the LRO state with the strongly 
reduced magnetic moment of about 0.4 fis- However, 
some of the previously reported experimental results re- 
quire further consideration. For example, more sensitive 
neutron-scattering experiments would be essential to de- 
tect or refute the LRO in CuNCN. Although the ^SR and 
NMR response might be characteristic of a partial inho- 
mogeneity, the origin of this effect could be purely mag- 
netic (e.g., the formation of a spin-glass state), purely 



structural, or both. We have shown that even stoichio- 
metric samples of CuNCN feature a sizable amount of 
defects that may affect the /iSR and NMR signals. An- 
other interesting result is the variable chemical compo- 
sition of CuNCN and the accommodation of foreign el- 
ements, such as oxygen and hydrogen, within the same 
structure type. Our study of the non-stoichiometric sam- 
ple 2 reveals the enhanced inhomogeneity of the magnetic 
state probed with /xSR and at least a partial destruc- 
tion of the LRO state. Therefore, intentional deviations 
from the ideal stoichiometry may be a feasible way to ap- 
proach the fully inhomogeneous magnetic grou nd state of 
CuNCN, similar to the one claimed in Ref. [l7| for the sto- 
ichiometric compound. 
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